Polaron Formation in the Three-Band Peierls-Hubbard Model for Cuprate 

Superconductors |cond-mat/9401025 



a^ 



> 

in 



o 



u 



J. Lorenzana 

Laboratory of Applied and Solid State Physics, Materials Science Centre, 

University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 



A. Dobrj 
Groupe de Physique Statistique, Universite de Cergy Pontoise, 47 Avenue des Genottes, 95806 Cergy Pontoise Cedex, France 

(February 1, 2008) 

Exact diagonalization calculations show a continuous transition from delocalized to small po- 
laron behavior as a function of intersite electron-lattice coupling. A transition, found previously at 
xj , Hartree-Fock level [Yonemitsu et al., Phys. Rev. Lett. 69, 965 (1992)], between a magnetic and a 

O^ ■ non magnetic state does not subsist when fluctuations are included. Local phonon modes become 

softer close to the polaron and by comparison with optical measurements of doped cuprates we 
conclude that they are close to the transition region between polaronic and non-polaronic behavior. 
^ ' The barrier to adiabatically move a hole vanishes in that region suggesting large mobilities. 
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The discovery of high To superconductors |]1| has trig- 
gered a renewed interest in strongly correlated systems. 
More recently the interplay between strongly correlated 
carriers and lattice effects has produced a lot of experi- 
^»\j mental and theoretical attention p] . Recent infrared (IR) 

^^ , absorption experiments show that added carriers produce 
substantial lattice relaxation around them showing up in 
phonon side bands and typical shake-up bands related to 
Cu-0 stretching modes p[-|5|. Similar effects were previ- 
ously found in photodoped^,p| experiments. 
d The nature of doping states coupled to the lattice has 

G . been studied in inhomogeneous Hartree-Fock (HF) with 
I ' random phase approximation (RPA) fluctuations calcu- 
^ , lations ||8|-[l0(|. This approach predict the expermentally 
iv ' observed doping induced phonon side bands as a con- 
^ , sequence of polaron formation, however the competition 
IL/ '■ between the self-trapping effects and the kinetic energy of 
• I— I , the hole is absent. The same problem arises in a recently 
^ ■ introduced LDA-t-U scheme ^. 

In this work we have performed exact diagonalization 
and HF with RPA fluctuation calculations of the three- 
band-Peierls-Hubbard model describing the Cu-0 planes. 
The exact calculations conflrm the formation of a local- 
singlet-polaron (LSP) for large electron-lattice coupling 
strength (A), analogous to the lattice polaron state found 
in HF and LDA+U. However the mean held picture is 
qualitative wrong for small A. Contrary to the HF and 
LDA-I-U approaches a minimum electron-lattice coupling 
strength (Ac) is necessary to produce the polaron in ac- 
cordance with the situation in uncorrelated systems p2[ . 
Furthermore by extending the previous HF results H] 
and comparing with exact diagonalization we show that 
a previously found transition between a magnetic and a 
non magnetic state on Cu is a feature of the approxi- 
mations involved and have no counterpart with the ex- 
act results. We show also that as the polaron is formed 



in hole doped materials local softening of the phonon 
modes occur. This explains the doping induced phonon 
side bands and by direct comparison with experiments 
we can estimate how far from Ac a real material is. Our 
result suggest that doped cuprates are close to the cross- 
over region where polaron formation occurs. We find that 
in that region the barrier for adiabatic motion of a hole 
vanishes suggesting relatively large mobilities. 

We consider a three-band-Peierls-Hubbard model 
&jlO[ in which for simplicity Cu atoms are kept fix and 
O ions are allowed to move with displacements Ui in the 
direction of the Cu-0 bond. Holes have an on-site inter- 
action Ud on Cu and f/p on O, a Cu-0 repulsion Upd, on- 
site energies on Cu Ed, and on O Ep, and Cu-0 hopping 
t (0-0 hopping is neglected here). When an O moves in 
the direction of one Cu with displacement \ui\ the corre- 
sponding Cu-0 hopping changes by a\ui\. Opposite signs 
apply when the O moves in the opposite direction. Each 
Cu-0 bond has a spring with force constant K |^-[ic|] . An 

2 

electron-lattice coupling constant is defined as A = ^. 
Parameters are Ep — Ed = 3,Ud — S,Up = 3, Upd — 
l,tpd = l,K = 32tpdA-'^ |Jl|l. Details of the method 
are presented in a preliminary report p3[ . 

The formation of the lattice polaron can be understood 
by repeating Zhang and Rice's |Q (ZR) construction in 
the case in which the O ions are allowed to relax adia- 
batically in a symmetric breathing mode along the Cu-0 
bonds with displacements \ui\ = u. Following ZR con- 
sider first a single CUO4 cluster with two holes. The 
energy of the local singlet state as a function of u is 
E{u) = e-lSu + AKpu"^, where e = A - 8t2(i + .^^-^) is 
the energy of the local singlet in absence of lattice relax- 
ation. P — 16at{^ + jj \^ ) is a renormalized coupling 

with the lattice. Kp = K — 2a^ (^ + jj-^zte. ) ^^ ^ renormal- 
ized spring- force-constant. The corresponding phonon 
frequency of the breathing mode uii = 2Kp/M , with M 
the O mass, is smaller than the one for the "undoped 
case" (only one hole) loq — 2Ko/M, with Kq — K — ^. 
A renormalized electron-lattice coupling strength can be 
defined as Xp = -^-^. Large \p implies strong tendency 
for polaron formation. The expression for (3 (note the 
large factor in front which comes from coherence effects) , 
shows that the strong covalency of ZR wave function fa- 
vors polaron formation and competes with A and Ud- 
We can find equilibrium positions for the lattice displace- 
ments vP — g^. The resulting state is a LSP showing 
local softening of the lattice breathing modes. Note that 
in the electron doped case (no holes) a hardening is ob- 
tain. We show below that IR active modes in the hole 
doped case behave in a similar manner in accordance with 
previous RPA results |8| . When the ZR singlet is allowed 
to hop, a competition will arise between the self-trapping 
and the delocalization effects. 

To further substantiate these findings we have per- 
formed a Lanczos diagonalization and HP plus RPA 



study of the cluster in the inset of Fig. |[ The O ion 
mass was taken to be infinity (adiabatic approximation), 
the O displacements where determined by iterating self- 
consistent equations for unconstrained O^s displacements. 
As we show below this is particularly important close to 
Ac. In Fig. 1^ we show the O displacements as a function of 
A in the two approaches. The higher branch corresponds 
to an O neighbor to Cui (see inset), so it represents the 
equilibrium displacement in the site of the polaron (Up) 
and the lower branch correspond to an O neighbor to 
Cu4 (w^). For A < Ac both branches coincide in the 
exact result and there is no polaron formation whereas 
for A > Ac the particle self-traps around Cui and the 
displacement there grows rapidly. The non zero lattice 
displacement for A < Ac indicates a tendency to charge 
density wave formation (CDW) in which the O's around 
Cui and Cu4 slightly contract and is discussed elsewhere 
|l3| . Contrary to the exact results the HF calculation 
shows \c,HF = 0. This is a consequence of the fact that 
in this approach the particle is already self-trapped with- 
out coupling to the lattice ||q- [l(]| , [l5|| forming a magnetic 
polaron state reminiscent of a local singlet |1J] . The same 
is true in LDA-I-U. Hence the competition between the 
self-trapping effects and the translational motion of the 
hole is absent in these mean field approaches. Note that 
for large A, when the symmetry is broken in both ap- 
proachs HF and the exact result are in good agreement. 
This is because the displacements are determined by the 
one body density matrix which is optimized in HF. In 
contrast, as we show below, two body correlations are 
badly reproduced. 

It was found in HF |P] that as A approachs a cross- 
over value A* a transition occurs to a non-magnetic state 
in Cui. This cross-over is signaled by an increase in the 
double occupancy of Cui towards the non-magnetic limit 

(("-TCui?HCui)noninag. = V^uJ^) ^ith TT-Cui = ("TCui + 

n|Cui) as shown in Fig. y by the triangles. In our case 
A* ^ .6. We have calculated the evolution with A of the 
exact double occupancy at this Cu. The result is shown 
by the open squares in Fig. ||. For A < Ac one should be 
aware that HF breaks translational symmetry in a dif- 
ferent way to that of the exact result, so it is difficult 
to compare both solutions in detail. For A = this can 
be partially solved by averaging over the four possible 
locations of the mean field state ||l^,|l6[ (large triangle in 
Fig. H). The result overestimates the double occupancy. 
For A > Ac both the exact and HF results break transla- 
tional invariance in the same way. It is clear in that limit 
that the exact double occupancy is well below the non- 
magnetic limit. We note also that the sudden increase 
at Ac does not refiect in the total Cu double occupancy 
that remains constant and it is rather an effect due to the 
rearrangement of the charges. To identify the source of 
the discrepancy we add fluctuations to the HF double oc- 
cupancy by calculating the RPA correlation energy and 
using the Hellmann-Feynman theorem. Remarkably we 



find that RPA partially restores translational invariance 
so that, contrary to the HF case, there is little difference 
for A '^ and X ^ Xc between the double occupancy in 
Cui alone and the spatially averaged value (see Fig. ||). 
The intermediate region is discussed below. In the pola- 
ronic region we show the result at Cui whereas in the de- 
localized region the spatial average is used (small circles). 
We also show the result at Cui for A = (big circle) . As 
A* is approached fluctuations become too big and the 
RPA breaks down. This is obvious in Fig. where at A* 
the RPA overestimate the correction and a non-physical 
zero value of the double occupancy is obtained. If it is 
calculated on Cui alone it becomes negative. In the pola- 
ronic region the correction is overestimated but it is clear 
from both approaches that the Cu remains magnetic. 

One should note that whereas HF shows a transition 
between a magnetic and a non-magnetic state in Cui the 
two solutions in HF plus RPA above and below A* rep- 
resent the same physical situation in two different ways. 
Since the effect of increasing A is to relax the O towards 
the Cu we can use A as a measure of the local cova- 
lency in the polaron. Consider the CUO4 cluster with 
two holes as before. In the limit in which covalency dom- 
inates (A > A* ) the HF wave function has the form 
1^^) = nj<t>pPi + <t^ddim whereat = i^^tiP]. and 
the sum runs over the four O. p- (djj.) creates a hole 
on O (Cu). This wave function has spin zero but large 
double occupancy on Cu given by 0^. For an occupancy 
on Cu of order one this correspond to ~ 1/4 as is the 
case in Fig. g for large A. RPA introduces correlations 
and this reduces double occupancy as is clear in Fig. 
for large A. The net result is a singlet with small dou- 
ble occupancy which looks very much like the ZR wave 
function. On the other extreme when covalency is not 
so dominant [A < A*] the HF wave function has broken 
symmetry and the single particle orbitals depend on spin, 
that is \HF) = Ylai^t'p.'rPi + 0d,adL)|O). This describes 
a magnetic state on Cu, say with up spin, and on O with 
opposite spin when 0^ ^ < 0d t and c/)^ ^ < (j)^ ^. The HF 
wave function is not a singlet now but RPA corrects for 
that through terms of the form Jcu~o^{S^S~ + SjSp) 
present in the residual interaction. Here Jcu-O is the 
Cu-0 exchange. RPA also reduces further the double 
occupancy and the net result is, as before, similar to a 
ZR wave function. Now is evident that the transition 
between these two HF solutions has no physical mean- 
ing but is a consequence of how HF plus RPA reproduce 
the physics in two different situations. An analogy can 
be made with the problem of a magnetic impurity in a 
metal. Anderson's classical treatment [|7| shows also a 
magnetic and a non-magnetic solution at HF level which 
can be associated with the magnetic and non-magnetic 
state of the Cu. Note the similarity between Fig. 5 in 
Ref. and Fig. 1(b) in Ref. g. Our result shows the 
danger of using mean field alone as a criteria to distin- 



guish between magnetic and non-magnetic states. 

To study how the LSP moves through the lattice we 
had parametrized the lattice displacement which interpo- 
lates between the LSP in Cui and in Cu4 |lj] . Note that 
the transition between nearest neighbors Cu is not al- 
lowed in a Neel ordered infinite lattice without producing 
excitations |jl9] . In Fig. a we show the total energy as a 
function of a collective coordinate (S) . The displacement 
are set as Ui = ^{u^ + u^) ± S and the signs are chosen 
in such a way that for A > Ac and S = ^{Up — u^) the 
configuration of the polaron sitting in Cui is obtained. 
For A < Ac there is a single minimum which corresponds 
to no polaron formation whereas for A > Ac there are 
two minima which correspond to the LSP siting in Cui 
or in Cu4. Note that at the transition point the barrier 
vanishes. Is essential for this result to consider a full re- 
laxation of all O in the cluster. If only the O around Cui 
are relaxed, the barrier never vanishes and the transition 
becomes first order like ||l3]. This is qualitatively differ- 
ent from a previous Holstein-Hubbard study Q . There 
coupling with O c-axis displacements was considered and 
the transition between the polaronic and non-polaronic 
state was found to be always first order like. 

Many experiments show local softening of IR modes 
upon hole doping ||3|,|6|,[7| or hardening for electron doping 
0. To study this effect we had made a frozen phonon 
calculation by moving in phase the two O at the right and 
left of Cui and computed the corresponding phonon fre- 
quency (wp) and compared with the result in the undoped 
case (loo) ( inset of Fig. ^ |2^]. The local softening effect 
is particularly well resolved in photodoped La2Cu04 [||J^] 
where the 708 cm~^ Cu-0 stretching mode gets bleached 
and a new band grows at 640 cm^^. This corresponds 
to wg/w^ ~ 1.22) and a A - 10% higher than Ac ^. 
This result, and the fact that effective masses in real 
materials are not very large, point to a situation where 
the real A in the cuprates is close to Ac. The physics 
close to that point is highly non-trivial for a finite value 
of the ion mass ll^p^ . The transition is not any more 
sharp but is washed out by quantum fluctuations and 
Ac becomes a cross over value. The renormalization of 
the effective mass of quasiparticles is moderate and not 
exponentially large like in the large A limit. An appeal- 
ing scenario is that for small doping polarons are formed 
with A close to Ac. This will result in a small activation 
energy for thermal diffusion of the hole as is observed in 
lightly doped cuprates y,^. For higher doping, polaron 
formation is inhibited due to phase space restrictions and 
a more complicated object occurs sharing characteristics 
of the self-trapped and non-self-trapped state and hence 
more metallic behavior. 

Very recently an infrared band at .1 eV in both hole 0] 
and electron hole [Q doped compounds has been identi- 
fied as arising from shake-up process of a polaron phonon 
cloud. In the electron doped case H] it was possible 



to separate the phonon contribution which appears as 
higher harmonics of the Cu-0 stretching modes consid- 
ered here, indicating the relevance of these modes in 
forming the polaron. Similar structures are present in 
the hole doped data p|,p^. 

We believe that these effects can be generic to a wide 
class of materials and provide a framework for classifi- 
cation. For example low mobilities in compounds like 
La2_a;Sra;Ni04+5 ]2^ , p5| , pl| . can be understood as aris- 
ing from a value of A >> Ac- In the former the self trap- 
ping probably inhibits the softening of the MIR band p^] 
observed in superconducting compounds ||27[] fo r which a 
physical picture was recently presented |l(j]l5t] . On the 
other hand we suggest that for La2-a;Sra:Cu04 probably 
A ~ Ac where interesting physics can arise. 

In summary, we have found a continuous transition be- 
tween a LSP polaron and a delocalized state in a three- 
band-Peierls-Hubbard model. We showed that there is 
no quenching of magnetic moments as found in previous 
approachs and by comparison with optical measurements 
we suggest that cuprates are close to the transition re- 
gion between polaronic and non-polaronic behavior. This 
leads to a rich phenomenology where lattice anharmonic- 
ities coexist with presumably high carrier mobilities. 
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FIG. 1. Equilibrium positions for one oxygen ion sur- 
rounding Cui (higer branches) and for another surrounding 
Cu4 (lower branches) versus A. The triangles (squares) are the 
HF (exact) result. The inset shows schematically the cluster 
studied in the presence of a polaron in Cui . 



FIG. 2. Double occupancy versus A. Full small (big) 
triangles are the HF result on Cui alone (averaged over the 
four O), open triangles are the unpolarized hmit on Cui, open 
squares are the exact result on Cui and full circles the HF plus 
RPA result averaged on the four Cu's (A < Ac) and on Cui 
alone (A > Ac). We also show the averaged result for A = 
(full big circle). 



FIG. 3. Total energy as a function of the collective coordi- 
nate 5 for A = .7 (diamonds), 1.05 (crosses) and 1.4 (circles). 
The inset shows the ratio uj'p/oJq versus A. 



